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ABSTRACT 

We use autocorrelation to investigate evolution in flow flelds inferred by apply- 
ing Fourier Local Correlation Tracking (FLCT) to a sequence of high-resolution 
(0.3 "), high-cadence (~ 2 min) line-of-sight magnetograms of NOAA active re- 
gion (AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar 
Optical Telescope (SOT) aboard the Hinode satellite over 12-13 December 2006. 
To baseline the timescales of flow evolution, we also autocorrelated the magne- 
tograms, at several spatial binnings, to characterize the lifetimes of active region 
magnetic structures versus spatial scale. Autocorrelation of flow maps can be 
used to optimize tracking parameters, to understand tracking algorithms' sus- 
ceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied 
include: time interval At between magnetogram pairs tracked, spatial binning 
applied to the magnetograms, and windowing parameter a used in FLCT. Flow 
structures vary over a range of spatial and temporal scales (including unresolved 
scales), so tracked flows represent a local average of the flow over a particular 
range of space and time. We deflne flow lifetime to be the flow decorrelation time, 
T. For At > r, tracking results represent the average velocity over one or more 
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flow lifetimes. We analyze lifetimes of flow components, divergences, and curls 
as functions of magnetic field strength and spatial scale. We find a significant 
trend of increasing lifetimes of flow components, divergences, and curls with field 
strength, consistent with Lorentz forces partially governing flows in the active 
photosphere, as well as strong trends of increasing flow lifetime and decreasing 
magnitudes with increases in both spatial scale and At. 



1. Introduction 



Estimates of photospheric velocities can be combined with the ideal Ohm's l aw to deter 



mine the fluxes of magnetic energy and helicity across the photosphere (e.g., IChad 12001 



Demoulin and Berge ] l2003l:ISchucljl2006h. These fluxes probably play important roles in driv- 



ing the corona to flare (IWelsch et al.ll2009f) and to launch coronal mass ejections (CMEs), as 
well as coronal heating (ITan et al.ll2007| ). Sever al techniques can be used to estimate photo- 



spheric fl ows, including tra cking methods fe.g.. iFisher and Welsch l2008l : ISchuckll2006[ ). ball- 



tracking (jPotts et al.ll2004l ) and feature tracking (IDeForest et al. 



2om. Techniques for es- 



timating photosph eric velocities, however, are imperfect ( iRieutord et al.ll200ll : IWelsch et al. 
20071 : ISchuckll2008l ). so efforts to improve estimation techniques are ongoing. 



Autocorrelation of image sequences has been used to quantify evolutionary timescales of 
structures in the solar atmosphere. In this approach, a measure of correlation (e.g., linear 
or rank-order correlation) is computed between pixel values at some initial time ti and the 
corresponding pixel values at a later time ti + AT, where we refer to AT as the lag time. 
As image features evolve wi th increasing AT, mea sures of correlation computed this way 
typically drop. For instance, iTritschler et al.l ( 120071 ) investigated interi ietwork and n e twork 



flne structure of the quiet Sun at disk center in the Ca II K line. I Welsch et al.l (|2009[ ) 
autocorrelated 96-minute-cadence, line-o f-sight (LOS) magne tograms of a few dozen active 
regions imaged by the MDI instrument ( IScherrer et al.lll995l ) in full-disk mode (with ~ 2" 
pixels) over lag times ranging from the data's nominal 96-minute cadence to several days. 
They found decorrelation times for the magnetic fleld, deflned as a drop in the autocorrelation 
to 1/e ~ 0.37 , on th e order of 72 hours, with considerable variation between active regions. 



Welsch et al.l (120091 ) also estimated horizontal flows in their active region sample by applying 
tracking methods to pairs of magnetograms separated in time by a tracking interval. At, 
equal to the data's 96-minute cadence. The flows were independ ently estimated by using 
two distinct methods, Fourier local correlation track ing (FLCT: iFisher and WelschI |2008| ) 
and the differential affine velocity estimator (DAVE; Schuck 20061 )^ Apodization windows 
of 8 and 9 pixels, respectively, were used, so motions on smaller scales were averaged away 
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over the apodization window. By autocorrelating the x- and y- components of the estimated 
flows, they found flow lifetimes of two to three 96-minute steps on this ~ 12 Mm spatial 
scale. 



The autocorrelation approach employed by lWelsch et al.l (|2009[ ) was useful for estimating the 
lifetime of magnetic field structures and flows on the particular spatial scales they studied. 
But what are the lifetimes of magnetic and flow structures on different spatial scales? How 
do flow lifetimes vary as a function of length scale? And what are the lifetimes of flow 
properties, such as vortical or converging motions? From the general properties of turbulent 
flows like those operating at the solar photosphere, we expect structure in flow fields to be 
present over a range of spatial scales, with a range of lifetimes. 

Further, little effort has been undertaken to understand how variations in tracking cadence 
At, noise, and choice apodization window size affect flow in : 



these quantities influence estimated flows? lAbramenko et al.l (120111 ) used feature tracking to 



'erence. How do variations in 



follow bright points over a range of At's, and characterized their motion as sup er-diffusive. 



The av erage speeds implied by their diffusive model decreased with increasing At. IChae et al. 



( 120001 ) applied an LCT algorithm to chromospheric H-a images, and noted that At must be 
long enough for displacements to be large enough to be detected by LCT, but not so long 
that displacements are larger than the apodization window. They also found th at choosing 



window sizes that were too small increased random velocities due to noise. Chae et a. 



( I2OO4J ) found slower average speeds for longer At's. More recently, IVerma and Denkerl (120111 ) 



conducted an investigation of the effects of parameter selection on flows inferred by applying 
local correlation tracking (LCT) to high-cadence, high-resolution (0.11 "pixel scale) G-band 
images of a field of view containing sunspot umbrae, pen umbrae, and quiet-S un regions, 
observed by the broadband filter irn ager (BFI) of the SOT (ITsuneta et al.ll2008l ) aboard the 
Hinode satellite ( iKosugi et al.l 120071 ). The evolutionary time scales and spatial scales of the 
underlying image structures motivated their selection of tracking intervals between 15 and 
200 s, and apodization windows of 600 - 2400 km. They found decreases in flow speeds both 
as the apodization window size increased, and as flow maps were averaged over longer time 
periods. Because photospheric magnetic flux is long-lived, tracking methods can be applied 
to magnetograms separated in t ime b y much larger tracking intervals — e.g., longer than an 
hour, as done by lWelsch et al.l (120091 ). 



Here, we address these and other topics by using FLCT to track a nearly 13-hour sequence 
of high-resolution (0.16" pixels, with a ~ 0.32" diffraction limit), rapi d cadence (^ 2 min) 



LOS magnetograms of NOAA AR 10930, recorded by the SOT /NFI (iTsuneta et al.ll2008 



Suematsu et al.l l2008l : llchimoto et al.l l2008l : IShimizu et al.l |2008|) instrument aboard Hinode 
over 2006 December 12 - 13. This active region was the source of an X-class flare and 
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CME during our tracking interval: the flar e began around 13 -Dec-2006 02:15UT and peaked 
around 02:40 UT in GOES 1-8 A X-rays jTan et al.l ( 2009 ) also investigated flows in this 
active region, focusing on the relationship of penumbral flows and evolution to the flare, and 
Schrijver et al.l (120081 ) conducted a comparative study of coronal magnetic fleld extrapola- 
tions derived from vector magnetograms of this active region. 

This paper's primary goals are: flrst, to investigate the effects of choices of tracking parame- 
ters on inferred flows; and second, to characterize the lifetimes and other properties of flows 
as functions of spatial scale and magnetic fleld. We begin by describing the magnetogram 
dataset in more detail (Q, then describe the methods we used to estimate flows (^, present 
our analysis of the lifetimes of flelds (§!]) and flows (^5]), and conclude with a brief summary 
of our results and their implications 



2. Hinode SOT/NFI Magnetograms 

We studied a sequence of Hinode SOT/NFI Fe I 6302 Ashuttered magnetograms of AR 
10930 with 0.16" pixels, created from the Stokes V/I ratio in Level data, recorded between 
12-Dec-2006 at 14:00 and 13-Dec-2006 at 02:58. The USAF/NOAA Solar Region Summary 
issued at 24:00 UT on 12-Dec-2006 hsted AR 10930 at S06W21, meaning it was near disk 
center during this interval. At this wavelength, the diffraction limit of the Solar Optical 
Telescope's 50 cm aperture is 1.22A/a ~ 0.32". We show initial and flnal magnetograms 
from the sequence in Figure [H At that stage of the Hinode mission, a bubble present within 
the NFI instrument degraded image quality in part of the fleld of view. In these frames, 
the affected pixels are in an area between from about 300 - 900 in x and above about 900 
in y. To account for the bubble, we exclude pixels above y = 848 in all our analyses of 
flow properties and lifetimes (although this region is retained in some images shown in this 
paper). 

To convert the measured Stokes / an d V signals i n to pix el-averaged flux densities, we applied 



the approximate calibration used by llsobe et al.l (120071 ) , 



= ~0?^ (Mxcm-2) , (1) 



where -Blos is the estimated LOS flux density, C/ and Cy are counts in t he / and V 



i mages 



respectively, and Cy/0.798 gives the circular polarization. As noted by llsobe et al.l ( 120071 ). 
this linear scaling is only approximate. In particular, in the dark cores of umbrae, the 
decrease in I introduces nonlinearity into the V/I ratio. This effect can produce spurious 
weak fleld regions in umbrae. A more complex, nonlinear calibration would be required to 
overcome this artifact. Since our primary goal is to investigate the time evolution of structure 
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in this magnetogram sequence, we have not pursued efforts to correct such artifacts in the 
V/I magnetograms. 

A detail ed spectroscopic ana lysis of fields in this acti ve region using Hinode SpectroPolarime- 
ter (SP: lTsuneta et al.ll2008l ) data was undertaken by lSchrijver et al.l (|2008[ ). who created one 
vector magnetogram with approximately 0.63 "pixels near the middle of our tracking interval, 
around 21:00 UT.Reprojection employed in their reduction procedure limits the ability to 
directly apply those results with our magnetogram data set, but, as described in Appendix 
[K\ we manually co-registered and resampled this SP vector magnetogram to map vertical 
magnetic field Bz and field strength |B| from this magnetogram onto the (4 x 4)-binned NFI 
magnetogram closest to the midpoint of the 45-minute SP scan (at 21:52 UT). 

The AR was at about S06W21 toward the end of our observations, and NFI's field of view 
(FOV) subtends about 10 heliocentric degrees at disk center. These facts imply cosines 
between the line of sight and the local vertical direction within the AR ranging from 0.88 to 
0.96, corresponding to variations in pixel scales from foreshortening of a few percent across 
the FOV. Over the earlier part of magnetogram sequence, upon which much of our analysis is 
focused, these variations are slightly smaller. Consequently, we take these modest distortions 
as acceptable, and opt not to reproject the magnetograms prior to tracking. Reprojection is, 
ho wever, certain l y app ropriate for tracking studies over much larger fields of view, e.g., that 
of IWelsch et al.l (120091 ) . who estimated velocities simultaneously across fields of view out to 
45° from disk center. 

Apart from three abnormally large time steps of approximately 600 s each, and two relatively 
small time steps of 26 s each, the 380 frames we analyze have a cadence of 121.4 ± 1.2 s. 

To study frame-to-frame changes in the magnetic field, we must first remove artificial evolu- 
tion due to whole-image shifts (rigid motion of each image) between frames, due primarily 
to variations in spacecraft pointing. To accomplish this, we first determined the frame-to- 
frame shift between the central regions of each pair of images (cropping the outer 10% of 
pixels along each edge), using FFTs to compute the frame-to-frame cross correlation func- 
tion and a second-order Taylor series appr oximation to estimate the displacement vector of 
the function's peak, to sub-pixel accuracy (IFisher and Welschll2008l ). We then used Fourier 
interpolation to estimate pixel values implied by a shift of the negative of the cumulative 
displacement from the initial frame in the sequence. Any wrapped pixels arising from this 
shift, due to the assumed periodicity, were then zeroed (as may be seen in a thin strip along 
the bottom of the right panel of Figured]). In Figure [2], we show the cumulative shifts in 
X and y over the duration of the magnetogram sequence. Large shifts in image alignment 
occur with an approxirn ate periodicity similar to that of Hinode's 98-minute orbital period. 
(IKatsukawa et al.l 120101 ). 
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Fig. 1. — Left: Initial Hinode/NFI LOS magnetogram (white is positive Hux, black is negative) from the 
13-hour data sequence we investigated, at full resolution (0.16 "pixels) with the saturation level set at ±500 
Mx cm^^ . Right: Final magnetogram in the sequence. 



In Figure [2l we also show the total unsigned flux and negative of the net flux (there is more 
negative LOS flux than positive; but the right plot axis is positive, so we changed the sign 
of the net flux), which both exhibit periodicities on a timescale similar to that of the whole- 
frame shifts. Analysis of Fourier spectra of the time series in pixels (not shown) shows a 
signiflcant excess power near the orbital frequency, but little evidence of helioseismic P-mode 
leakage into the magnetogram signal. Given the large Doppler shift from orbital motion, 
observed periodicities could arise from leakage of Doppler shifts induced by spacecraft orbital 
motion into the Stokes V measurement. Thermal variations in the instrument might also be 
responsible. 

As discussed in more detail below, local fluctuations in fleld strength from noise can mimic 
changes from true magnetic evolution, implying estimated flows will be inaccurate in regions 
where the measured flelds ar e near the magri e tograp h noise level. We estimated this noise 



level using the procedure of iHagenaar et al.l (119991 ). First, we spatially binned the data 



2 X 2 to 0.32"-wide pixels, to approximate the SOT diffraction limit of 0.3". Next, we 
created histograms of pixel values in each frame in the sequence, in bins of 3 Mx cm~^. 
Then, assuming the core of the distribution in signals arises purely from noise, we flt a 
Gaussian, exp[—{B — i?o)^/(2a^), to the central region of the histogram, ±10 Mx cm~^. 
Over the course of the tracking interval, fltted values of the width parameter a for single 
magnetograms ranged from a high near 17 Mx cm~^ at the start of the run to lows near 
7 Mx cm~^ toward the end of the run. For simplicity, we therefore take 15 Mx cm~^ as a 
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Fig. 2. — Left axis: Cumulative frame-to-frame shifts (in pixels) in the x (red) and y (blue) directions be- 
tween pairs of magnetograms in the NFI sequence that we analyzed, with respect to the initial magnetogram. 
Right axis: total unsigned flux (dashed) and the negative of net flux (solid). Variations in flux probably 
arise from crosstalk due to Doppler shifts induced by spacecraft orbital motion into the Stokes V signal. 



uniform estimate of the noise level. 



3. Estimating Flows 



Assuming magnetic diffusi yity is negligibly small over t he time interval At between a pair of 
successive magnetograms, iDemoulin and Bergerl (120031 ) argued that the apparent motion of 
magnetic flux from one magnetogram to the next over At could be described by a "footpoint 
velocity" u, that is related to the plasma velocity v by 



VnBh 



(2) 



where n denotes the unit vector normal to the surface, the n subscript denotes the normal 
component of a vector, and the h subscript denotes the horizontal components of a vector, 
which are perpendicular to fi. This assumption can be used to express the normal component 
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of the magnetic induction equation in terms of a continuity-like, finite-difference equation 

AB 



At 



■ (uBn) = 



(3) 



where ABn = Bn{tf) — Bn{ti) is the change in the normal magnetic field between the initial 
time ti and the final time tf, and is approximated with finite differences in the spatial 
coordinates orthogonal to the normal direction. The analogy with the continuity equation 
is not precise, because unlike number or mass densities, which are always positive, magnetic 
flux can be positive or negative. So, unlike total particle number or total mass in a volume, 
neither total signed nor total unsigned flux is conserved in a magnetogram: existing opposite 
polarity fluxes can cancel, or new bipolar flux can emerge. 

Demoulin and Berger argued t hat t racking r nethqds app lied to magnetograms esti- 

mate u, instead of (cf., IChaell200l[ ). but ISchuckl (120051 l2008l ) argued that some tracking 
methods return a biased estimate, u, with the contribution from v/j weighted more heavily 
than in the true u. The term "footpoin t velocit y " imp licitly refers to extensions of magnetic 
flux above the magnetogram surface. IWelschI (120061 ) used the alternative "flux transport 
velocity" to refer to u without reference to field structure above or below the magnetogram. 
For purposes of tracking, we assume the LOS field, Bi, derived from the NFI observations, 
approximates the normal field Bn- 



We used the FLCT code (jWelsch et al.ll2004j : iFisher and Welschll2008[ ) to estimate flux trans- 
port velocities for this dataset. The FLCT code uses Fourier cross-correlation of Bn{ti) with 
Bn{tf) near a given pixel, say one centered at {xj,yk), to determine the two-component dis- 
placement AR of magnetic field structure between the magnetograms near [xj, Uk)'-, division 
by At then results in a velocity estimate at {xj,yk)- The correlation is localized by first 
weighting each magnetogram by a Gaussian windowing function, exp(— r|;,/cr^), where rjk is 
distance from {xj,yk), and a is a free parameter. Optionally, low-pass filtering, with a Gaus- 
sian roll-off, can be applied to the Fourier spectra p rior to performing the cro ss-correlation. 
In tests applying shifts to quiet-Sun magnetograms, iFisher and WelschI (120081 ) found roll-off 
values of 0.2A;max - O.Sfcmax improved recovery of the applied shifts, where fcmax is the highest 
wavenumber each dimension. In this work, we used a roll-off wavenumber of 0.5 fcmax, which 
corresponds to relatively weak Fourier filtering. 

Like FLCT, other tracking methods typically apply an apodization or windowing function, 
the spatial extent of which is a free parameter, to the initial an d final magne tograms, and 
find a velocity that extremizes some functional over this window ( ISchuckl |2006[ ). Given data 
with sufficient cadence, the interval At between pairs of magnetograms to be tracked can also 
be varied. Tracking methods therefore estimate velocities from the effects of actual velocities 
averaged both spatially over the windowed region and temporally over the time interval At 
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between magnetograms. This averaging implies that flows on spatial scales smaller than the 
windowing function, and with lifetimes shorter than At, cannot be recovered in the flow 
estimation process. 

Beyond these algorithmic aspects of the tracking process, it can be also seen that noise 
and other artifacts in magnetogram measurements will contribute to the difference AB^ in 
equation ([3]), implying such effects can influence estimated velocities. 

These considerations motivate efforts to understand the interplay between windowing and 
tracking cadence and flow spatial scale and lifetime, as well as the effects of magnetogram 
noise. To investigate these relationships, we estimated flux transport velocities from pairs of 
magnetograms that were separated in time by eight varying intervals, with At G {2, 4, 8, 16, 
32, 64, 128, 256} minutes, that were both unbinned and binned to six varying macropixel 
sizes Ax G {2, 4, 8, 16, 32, 64 } in terms of the data's original 0.16" pixels, with four 
choices of apodization parameter a G {2, 4, 8, 16 }. Pixel values were simply averaged 
when rebinned, which does not exactly mimic how magnetographs with either worse spatial 
resolution or poor weak-field sensitivity would detect magnetic flux. 

Since we are interested in characterizing the effects of noise on estimated flows, we tracked 
all pixels above a threshold of 15 Mx cm~^, which is both near our estimate of the noise 
level, and includes the spurious weak-field regions in the sunspot cores. 

4. Lifetimes of Magnetic Structures 

To provide context for the time evolution of flows on varying spatial scales, we autocorrelated 
the magnetograms themselves, at several spatial resolutions: with 0.16" pixels, and then 
with data rebinned to macropixels Ax G {2, 4, 8, 16, 32, 64} pixels on a side. While each 
magnetogram is a 2D array, we computed the correlations as if both were ID arrays. All 
pixels outside of the bubble region were used, not just those above the noise threshold used 
for tracking. The colored curves in Figure [3] show rank-order correlation coefficients for 
the autocorrelations at each spatial resolution for increasing lag time. In addition, the black 
dashed and solid curves show frame-to-frame linear and rank-order correlations, respectively, 
for the 2x2 binned data. Note the limited range of the y-axis. 

The high frame-to-frame correlations imply that the magnetic field does not evolve much 
from one magnetogram to the next. It can also be seen that frame-to- frame linear correla- 
tions hover around unity, but frame-to-frame rank-order correlations are significantly lower, 
implying rank-order correlation is a more sensitive test of differences between magnetograms. 
This plot also shows that the lifetime of this active region's magnetic structure, over the en- 
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Fig. 3. — Colored curves show rank-order autocorrelation coefficients with increasing lag times, for magne- 
tograms at several spatial binnings: red is for 0.16" -pixel data; orange, yellow, green, aqua, blue, and purple 
are for data rebinned 2 x 2, 4 x 4, 8 x 8, 16 x 16, 32 x 32, and 64 x 64, respectively. The black dashed and 
solid lines show frame-to-frame linear and rank-order correlation coefficients, respectively, for data binned 
2x2. Note the limited range of the vertical axis. Downward spikes before 500 minutes arise from frames 
with many cosmic ray hits. The effects of the X-class flare can be seen as small blips near 743 minutes. 

tire magnetogram field of view, is mucli longer tlian tlie duration of the dataset tliat we 
analyze. 

Several artifacts are also present in these time series. Downward spikes in correlations 
correspond to frames in which isolated pixels contained large fluctuations, probably due to 
hits by cosmic rays. Despiking the magnetograms can remove these spikes, but here we are 
interested in the effects of such artifacts on autocorrelations. Larger binnings average out 
these effects. In terms of lag times plotted here, the start time of X-class flare that began 
around 02:28UT on 2006-Dec-13 corresponds to 743 minutes, and its effects can be seen as 
small wiggles in some autocorrelation curves near the right end of the plot. 

Distinct types of magnetic regions — notably, sunspots, plage, and quiet Sun — all influence 
the autocorrelations in Figure |3l By investigating the lifetimes of magnetic structures on 
smaller spatial scales, we can better understand the lifetimes of magnetic structures in each 
type of region. Accordingly, in Figure HI we show autocorrelations (red curves) for -Blos 
in (32 X 32)-binned subregions of (2 x 2)-binned data (so each subregion is (64 x 64) full- 
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resolution pixels). In the same figure, we also show one-parameter fits (blue curves) to the 
decorrelation in each subregion, assuming the autocorrelation function decays exponentially 
with lag time AT, exp(— AT/r), with the decay constant r as the only free parameter. In 
each subregion, we take the fit decay constant, r, to be the lifetime of the structure in that 
subregion. We define the occupancy of a subregion to be the percentage of pixels above our 
15 G noise threshold. Only pixels with a median occupancy, over the dataset's duration, of 
at least 20% were fit. Background gray contours in the figure show 50 G and 200 G levels 
of I -Blob I from the initial magnetogram (in full- resolution pixels), so variations in lifetimes 
can be understood in the context of active region structure. 

We then investigated the dependence of magnetic lifetime on field strength. To start, we 
computed a measure of subregion field strength by boxcar averaging flux densities from 
I-BlosI from the five initial magnetograms in our sequence at (2 x 2) resolution, then taking 
the absolute value, then computing the spatial average of this unsigned, boxcar-averaged 
flux density over each subregion. We chose to use this measure of initial field strength to 
address the question: "For a given current field strength, how long will a structure persist?" 
The relatively short fitted lifetimes of some magnetic regions suggests that averaging over 
longer intervals — e.g., the length of the entire run — would imply averaging over several 
magnetic turnover times. It can be seen in Figure H] that some of the fits represent the 
evolution poorly, but many fits match the observations well, given the single free parameter. 
Of the 195 subregions with at least 20% above-threshold occupancy, only four have fitted 
lifetimes that are less than the two- minute nominal cadence (see, e.g., the two subregions 
to the left of the top-right subregion). We ignore these too-short fits as pathological. In 
Figure [5], we plot the fitted lifetimes as a function of these spatially and temporally averaged 
subregion fields (blue +'s). A linear fit between the logarithms of each returns a power- 
law exponent of 0.66. The solid blue line shows the fit, while the dashed blue lines reflect 
uncertainty in the fitted slope, based upon computed from the fit. This trend is consistent 
with the idea that convection disr upts magnetic field structures, but stronger magnetic field s 



more strongly inhibit convection (ITitle et al.l Il992l : iBerger et al.l Il998l : IWelsch et al.l |2009| ) 



so structures persist longer in stronger-field regions. Despite the trend of increasing hfetime 
with field strength, we note some high-field-strength regions have short lifetimes. Relatively 
rapid evolution in these high-field strength regions might play some role in flare activity, but 
pursuing this question is beyond the scope of this paper. 

To investigate the dependence of this field-strength-versus-lifetime relationship on spatial 
scale, we repeated this procedure with the subregion size reduced to 16^ (2 x 2)-binned 
pixels (autocorrelations and fits at this resolution are not shown). Lifetimes as a function of 
subregion- averaged field strength are plotted in red in Figure [5l excluding 28 subregions with 
fitted lifetimes < 2 min (of 733 subregions with at least 20% above-threshold occupancy). 
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A linear fit between the logarithms of field strengths and lifetimes returns a power-law 
exponent of 0.69 in this case. The solid red line shows the fit, while the dashed lines reflect 
uncertainty in the fitted slope. The trend of lifetime versus field strength did not change 
much for subregions smaller by a factor of two. 

Recalling that our estimates of |-Blos| suffer from non-monotonicity in the Stokes V/I ratio, 
we also characterized the lifetimes of NFI structures as functions of |-Bz| and |B| from the SP 
data, in subregions averaged to the same spatial resolution (16^ regions of 2 x 2 pixels). Note, 
however, that the SP data were recorded near the middle of our observing run (between 20:30 
and 21:14 UT), while the magnetic lifetimes were computed from the beginning of our run. 
Also, our treatment takes no account of spatial filling factors of magnetic structures within 
resolution elements. Despite these uncertainties in the SP field strengths, fits of lifetimes as 
a function of (not shown) were similar to those for |-Blos|- The solid green line in Figure 
[5] shows the fitted slope (dashed lines reflect uncertainties in the fitted slope) of lifetimes 
as a function of |B|. The slightly weaker dependence of lifetime on |B| compared to |-Blos| 
might arise from either: vertical magnetic fields inhibiting convection more efficiently than 
horizontal magnetic fields; or rapid magnetic evolution near polarity inversion lines (PILs) 
of the normal field, regions where the magnetic field is predominantly horizontal and where 
flux emerges or submerges; or both. We note that the range of |B| is greater than for either 
of the other tw o quantities, and this l onger run will naturally decrease the fitted slope for 



the same rise. iDemoulin and Bergerl (l2003l ) suggested that in regions with predominantly 



horizontal magnetic fields, apparent motions of vertical magnetic flux could be rapid. 

We remark that autocorrelation analyses from observed magnetograms like those presented 
here could be compared to autocorrelat ions from synt hetic magnetograms extracted from 



simulations of magnetoconvection (e.g., iRempell l201ll ). as a test of statistical consistency 
between the observations and simulations. 



5. Flow Lifetimes 

In this section, we first consider selection of tracking parameters, and then investigate life- 
times of estimated flows versus spatial scale and magnetic field strength. 



5.1. Selection of Tracking Parameters 



Welsch et al.l (l2004f) n oted that flows derived by tracking alone only approximately solve 
equation ([3]). ISchuckl (120061 ) suggested that, in the presence of magnetogram noise, exactly 
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solving the induction equation is probably not optimal. IWelsch et al.l (120071 ) suggested that 
for fixed At (i.e., imposed by available data), a could be chosen to optimize agreement 
between the observed ABn/At and the estimated Vh ■ (u-B„). When observations are fre- 
quent enough that At can also be varied, considerations beyond optimizing consistency with 
equation ([3]) are useful. 



In addition, IWelsch et al.l ( 120111 ) demonstrated that enforcing strict consistency with equa- 
tion ([3]) in deriving fiows is probab ly unwise. Using 96 -minute cadence, full-disk magne- 
togram sequences observed by MDI (jScherrer et al.lll995l ). they characterized magnetic evo- 
lution in three active regions near disk center. In addition to using FLCT to estimate flows, 
they also numerically derived ideal electric fields to exactly solve Faraday's law, from a 
Poisson equation, with a source given by their estimate of t he radial magnetic field evo- 
lution [A-Bff/At], f or the inductive part of the electric field (IWelsch et al.l 120071 ). Where 
Welsch et al.l (120091 ) found fiow fields derived by either FLCT or DAVE from 96-minute ca- 
dence magnetog rams to be co r relate d from one 96-minute interval to the next, the "Poisson 
fiow fields" that lWelsch et al.l (120111 ) derived to be exactly consistent with ABr/ At decorre- 
lated from one frame to the next. Assuming random fluctuations (due to, e.g., noise present 
in the estimated flows), it can be shown that, if evolution in the flow is negligible (so the 
only difference between the flow maps is due to random errors in the estimated flows), then 
even a relatively high correlation coefficient (e.g., 0.9) implies signiflcant fluctuations present 
(respectively, 33% of the variance in the flow component). It can also be shown that even 
allowing for evolution in the flow only moderately reduces the fraction of decorrelation that 
must be ascribed to random fluctuations in the flow estimates. The fact that FLCT and 
DAVE flows did not decorrelate between frames, but the Poisson flows did, implies the latter 
contained signiflcant random fluctuations, and were therefore probably strongly influenced 
by magnetogram noise. This suggests that using equation ([3]) as the sole guide to determining 
tracking parameters has limitations. 

We therefore consider the effect of the time interval At between each pair of tracked mag- 
netograms on the relationship between flux transport velocities u from tracking and plasma 
velocities v in equation ([2]). Essentially, if At is either too short or too long, then equation 
([2]) is invalid. To explain why, we consider each case separately. 

When the time difference At between magnetograms is too short, very little magnetic evolu- 
tion has occurred (as in the black curves in Fig. [3]), so most of the difference Ai?„ between 
initial and flnal magnetograms that appears in equation ([3]) is due only to noise. Conse- 
quently, flux transport velocities inferred from ABn/At will be spurious. For instance, a 
transient spike in the fleld strength of a single pixel in a weak unipolar region would lead 
to inference of a high-speed converging flow surrounding that pixel, consistent with the sud- 
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den field increase, followed by a diverging flow, consistent with the field decrease after the 
spike disappears. In general, these noise-induced flows will be unrelated to actual plasma 
velocities. Consequently, if At is too short, then noise is a large component of magnetogram 
difference ABn, and we refer to flows estimated in this regime as noise-dominated. 

When the time difference At between magnetograms is too long, displacements of magnetic 
flux might exceed a pixel length, meaning the finite difference approximation used in equation 
(E]) could be inaccurate. In this regime, which we refer to as displacement-dominated, flux 
transport velocities estimated by tracking are more properly described as time averages of 
true flux transport velocities, which acted and evolved over timescales shorter than At. It 
should be noted that, although estimated flows in the displacement-dominated regime do 
not correspond directly to plasma velocities, they still have a clear physical meaning: they 
represent the time-averaged flow over timescales of a flow lifetime or longer. We expect that 
motions derived for At longer than a flow lifetime will underestimate true flow speeds: the 
true path of a feature executing zig-zag motion can be resolved with short enough At, but 
if At is too long to resolve each step on the path, only the (shorter) net displacement is 
recovered (cf., the full path length traversed over that time interval). Hence, a given At can 
average over the effects of shorter-lived flows. 

It would be convenient if an ideal choice of At could be made a priori. As previously noted, 
however, flows with a range of speeds, and time and length scales are present at the solar 
photosphere (e.g., granular flows with speeds ~ 1 km sec~^ with lifetimes of a few minutes 
and ~ 1 Mm length scales, and supergranular flows with speeds a few times 10^ m sec~^ with 
lifetimes of about a day, and length scales of ~ 10 Mm). Poor choice of At can, however, 
be determined by post facto analysis of tracking results, meaning trial runs can be used to 
identify suitable timescales. 

We flrst consider the noise- dominated regime, in which correlations between each velocity 
component from successive flow estimates will be small. Frame-to-frame correlations can 
be improved by increasing At = tf — ti, and by averaging successive magnetograms around 
ti and tf to synthesize initial and flnal images. Figure [6] shows frame-to- frame, rank-order 
correlations for the x- and ?/-components of estimated velocities, with a flxed at 4 pixels but 
for several At's, using either unaveraged (dashed) and averaged (solid) magnetograms for 
B{ti) and B{tf). Low correlations for short values of At imply that noise has degraded the 
velocity estimates. For At = 2 min, the flow flelds are even weakly anticorrelated. Increas- 
ing At clearly increases frame-to-frame correlations, because real magnetic fleld evolution 
becomes more signiflcant than spurious evolution from noise. Saturation occurs, however: 
doubling At from 8 to 16 minutes did not greatly improve frame-to-frame correlations. As At 
increases further, into the displacement-dominated regime, frame-to-frame correlations can 
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approach unity, since displacements averaged over very long time intervals do not change 
much if the intervals' endpoints are shifted slightly. To average the magnetograms, a 3- 
magnetogram boxcar was used for At = 2 and 4 min, while a 5-magnetogram boxcar was 
used for At = 8 and 16 min. For cadences of 2, 4, and 8 minutes, there is overlap in data 
used for B{ti) and B{tf). This probably leads to blurring of features from their moti on over 



the s umming interval, and could also lead to correlated noise in B{ti) and B{tf) (jSchuck 



20061 ). The improvement in frame-to- frame correlations from averaging magnetograms is 
substantial. Linear correlation coefficients (not shown) were both lower and less consistent 
across frames. We therefore employ rank-order correlations to compute all subsequent auto- 
correlations. We provide concrete examples of these effects in Figures [7] and |8l In the former, 
successive flow fields computed with At = 2 min show very little consistency between fiow 
fields. In the latter, successive fiows computed with At = 8 min are much more consistent. 
Generally, frame-to-frame correlations are higher with larger values of a (not shown). This is 
consistent with FLCT detecting fiows on larger spatial scales with larger a, and larger-scale 
fiows persisting longer in time. We discuss this issue in greater detail below. 

In the displacement-dominated regime, correlations between each velocity component from 
fiow fields separated by lag times AT much shorter than At can be high, but such correlations 
will be very weak when AT is on the order of At or longer. This implies substantial fiow 
evolution has occurred over the At between magnetograms. Phrased another way. At is 
longer than the estimated fiows' lifetimes. For estimated fiows to refiect actual plasma 
velocities, however. At cannot exceed the fiow lifetime. Displacement-dominated fiows are 
still physically meaningful, however: they quantify the average effects of fiows over the course 
of more than one fiow lifetime. 

Having demonstrated the advantages of averaging, we use fiows derived from magnetograms 
averaged with a five-step boxcar in our analyses of fiow lifetimes below. For convenience, 
we will sometimes refer to fiows inferred by FLCT as velocities, but the reader should keep 
in mind that FLCT's fiows should not necessarily be interpreted as estimates of plasma 
velocities at the center of the time interval At. 



5.2. Estimating Flow Lifetimes 

We define the lifetime r of a fiow as the lag time AT at which a one-parameter fit to the 
autocorrelation drops to 1/e ~ 0.37. In Figure [H we plot autocorrelations for fiow estimates 
made with two choices of a and several choices of At. For At = 64 minutes, the lifetimes of 
the X- and y- fiow components are around 35 minutes for a = 8, implying that fiow estimates 
made with At = 64 minutes are displacement-dominated. In contrast, the lifetimes of cr = 8 
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flows estimated with At = 4, 8, and 16 minutes are all at least as long as At. Lifetimes of 
flows estimated with cr = 4 (dashed lines) are about half as long as cr = 8 estimates for At = 
8, 16, and 64 minutes. 

The autocorrelations in Figure [9] also suggest that flows estimated on larger spatial scales 
— using a larger value of a — have longer lifetimes. In fact, this is typical. In addition, 
longer flow lifetimes are also found when data are rebinned into larger macropixels, though 
there is a point at which rebinning does not increase decorrelation times. The top panel of 
Figure [10] demonstrates that choosing larger values of a or rebinning into larger macropixels 
are practically equivalent. The bottom panel of Figure [10] implies that, when flows are 
displacement-dominated (i.e., they decorrelate at lag times AT < At), longer-lived flows 
can be identified in rebinned data. We note that increasing the windowing parameter a can 
increase compute time for FLCT (and probably other tracking methods), while rebinning 
(holding a fixed) can decrease compute time, since fewer velocity estimates must be made. 

In Figure [11], we show autocorrelations of Ux and Uy in (32 x 32)-pixel subregions of mag- 
netograms with 0.32 "pixels (so binned 2x2). These flows were derived with At = 64 
minutes and a = 4 pixels. Background gray contours in the figure show 50 G and 200 G 
levels of I -Blob I from the initial magnetogram (in full-resolution pixels), so variations in flow 
autocorrelations can be understood in the context of active region structure. 

In weak-field regions, few pixels are above the threshold for tracking, so very few velocities 
are estimated. Consequently, the autocorrelations are essentially random. We therefore only 
fitted the autocorrelations in high-occupancy subregions, which we define to be those with a 
median (over the dataset's duration) occupancy of pixels above the tracking threshold of at 
least 20%, which gives large enough samples that the autocorrelations are stable. The fits to 
the magnetic field autocorrelation curves in Figure H] convey how well the single-parameter 
exponential fits match the autocorrelations for a range of lifetimes and autocorrelation pro- 
files. 

In Figure [121 we show fitted flow lifetimes {ux = red -|-'s and Uy = blue x's) for high- 
occupancy subregions of Figure [TT] as a function of the spatial average of the unsigned, 
boxcar-averaged LOS magnetic flux density in each subregion. Red A's and green V's are 
lifetimes of Ux components of the flow as functions of spatial averages of and |B|, respec- 
tively, and red dashed and green solid lines are fits to each. Most flow lifetimes are between 
about 7 and 30 minutes, confirming the expectation based upon the red, dashed curves in 
Figure [H] that these flows are displacement-dominated. There are statistically significant cor- 
relations both between flow lifetime and field strength, and between flow lifetime and median 
occupancy of above-threshold pixels. The rank-order correlation of field strength |B| with 
Ux (and Uy) flow lifetime, however, is a statistically significant 0.52 (and 0.51) for the 35 of 
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137 non-bubble subregions of SP data with occupancies of 1000 or more (of a possible 1024), 
implying the lifetime versus field-strength trend is independent of any occupancy-lifetime 
correlation. Linear fits to logarithms of lifetime as a function of field strengths return weak 
power-law exponents, below ~ 0.2 for each flow component. For flows with At = 8 min- 
utes and cr =4 pixels, correlations of lifetime with fleld strength are quite weak. For flows 
estimated with At = 256 minutes and a =16 pixels, fltted power-law exponents between life- 
time and fleld strengths {I-BlosI; \Bz\, |B|} range from 0.2 0.3, but fleld strength-lifetime 

correlations are weaker for Ux than for Uy. The east-west orientation of the PIL, combined 
with apparent rapid motion of vertical magnetic flux along the PIL (jPemoulin and Berger 
20031 ). could produce more rapid evolution of east- west flow components than north-south 
flow components. 

Correlations between flow lifetime in a subregion and magnetic fleld strength in that subre- 
gion imply that Lorentz forces drive a component of photospheric evolution. To the extent 
that convection is the primary driver of photospheric evolution, and magnetic flelds inhibit 
convection, there is no obvious reason why flow lifetimes should be longer in stronger-fleld 
regions. But because magnetic structures persist for longer times cales than flo ws, it is likely 
that forc es acting on those fle l ds, su ch as magnetic buoyancy (jParkerl Il955l ) or magnetic 
torques (ILongcope and WelschI l2000l ) are also longer-lived. Consequently, it is physically 
reasonable that in stronger-fleld regions, magnetic evolution is partially driven by longer- 
acting forces than the pressure gradients arising from convection. 

How does flow lifetime r depend upon flow spatial scale? A crude estimate of the turnover 
time of a convective eddy — expected to be related t o flow lifetime — is the eddy length scale 



L, divided by a characteristic speed uq. As noted by lBerger et al.l ( 119981 ) . however, "the term 
'characteristic speed' must be used only in relation to a given spatiotemporal resolution." So 
as a preliminary step in characterizing flow lifetime versus spatial scale, we flrst investigate 
flow speed versus spatial scale. What do our results reveal about typical flow speeds over 
varying spatial scales? We deflne a flow's spatial scale as the product of macropixel bin 
size and the windowing parameter a. As shown in Figure [TOl changes in both parameters 
seem to have similar effects upon flow lifetimes. Because we estimated flows for several 
macropixel sizes and choices of cr, we therefore have multiple sequences of flow estimates 
at a few spatial resolutions for each cadence At. We computed the average flow speed at 
each scale over all pixels and time steps for which velocities were estimated. In Figure [131 
we show average velocity as a function of spatial scale; the plot symbols are log2(bin size), 
in units of the data's native 0.16 "pixels, to better illustrate changes in average velocity as 
a function of bin size at a given spatial scale. (Values of curls and divergences of the flow 
pattern are also shown in this flgure, and are discussed below.) One basic property is evident : 
average speeds are slower for longer tracking intervals At, as was noted by lChae et al.l (120041 ) 
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and IVerma and Denkeii (120111 ) for both relatively long and short At, respectively. (This is 
discussed further below.) For At shorter than 16 min, a verage speeds peak at the small end of 
the range of spatial scale s, as found in previous work ( Berger et al.lll998l : IChae et al.l 12004 : 
Verma and Denkerl 1201 ll ). For longer At, however, average speeds peak at higher spatial 
scales: speeds at small-scales for large At are low. The large At implies that small-scale 
structures at tj and tf have evolved substantially. While it is plausible that FLCT would 
infer random displacements in this case, and that averaging over these random displacements 
(all positive) would produce a large mean speed, we do not observe this. Apparently, for a 
given At, the FLCT code is primarily sensitive to flows on a particular spatial scale. Also, for 
shorter At, average speed at a given spatial resolution depends strongly upon spatial binning, 
while for longer At, average speeds seem insensitive to spatial binning. This probably reflects 
different effects from the combination of actual evolution and noise on the apparent magnetic 
evolution convolved with spatial binning and apodization by a. For instance, large values 
of a presumably average over the effects of actual evolution on sub-cr scales. If so, speeds 
should be higher with smaller a at a given spatial resolution, as we find. We remark that the 
combined effects of cadence with binning and apodization preclude straig htforward testing 
our FLCT velocities for consistency with the Kolmogorov scaling (e.g., iThompsoru |2006| . 
p.68). 

We note that tests of FLCT with synthetic magnetograms from MHD simulations, in which 
the true vel ocities were known , showed that FLCT systematically underestimates speeds, 
by ~ 25% ( IWelsch et al.l 120071 ). so the absolute values of the average speeds we find here 
are probably systematically lower than estimates that would be found using other tracking 
codes. Further, because our velocities are derived from magnetic evolution, and magnetic 
flux both tends to inhibit convection and to concentrate in stable downflow regions between 
convective cells, it is reasonable to expect our average speeds to be systematically lower than 
tracking results derived from intensity images. We note that for A t of 2 or 4 min, our average 
speeds at the smallest spatial scale are close to those repor ted bv IVerma and Denkerl (120111 ) 
(see their Table 3). The speeds of FLCT flows derived by IWelsch et al. (l2009h . with At of 



96 min and a spatial scale of ~12 Mm were predominantly in the 50 - 100 m sec ^ range 
(see their Fig. 6), comparable to values we find for similar At and spatial scales here. 



In Figure [TU we show the average speed versus tracking cadence At, averaged over both 
space and time in all flow maps at a given spatial resolution and all spatial resolutions 
at each cadence. Effectively, this amounts to computing the average of speeds at each 
At in Figure [T3] — fitting the speeds in each panel to a horizontal line. Error bars in 
Figure [H] reflect the standard deviation in speeds in each panel of Figure [121 A linear fit 
to the logarithms of cade nces and average speeds returns an power-law exponent of -0.34. 
Verma and Denkerl ( 1201 ll ) computed average speeds for a range of much shorter At's, and 
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over the eight-fold increase in At from 60 to 480 sec, their average speed fell by about a 
factor of two. This implies a scaling near -1/3, similar to ours. (But they found an increase 
in average speed between Afs of 15 and 60 sec, and similar average speeds at 60 and 90 
sec.) They suggest decreasing speeds with increasing At might arise from the short lifetimes 
of features in the images they track. However, we find the same effect when tracking much 
longer-lived magnetogram structures, implying the limited lifetimes of features is probably 
not solely responsib l e for decr easing speeds with increasing At. We also note that both 
Hagenaar and Shind (120051 ) and lVerma et al.l (120111 ) found that individual, successive moving 
magnetic features (MMFs) followed preferential paths as they moved outward from sunspots, 
suggesting short-lived features can reveal longer-lived flow patterns. As noted above, a longer 
At can average over the effects of short-lived flows: displacements inferred with a longer At 
might not recover zig-zag motions that can be resolved with shorter At, with the shorter 
apparent path length at longer At implying a lower estimate of speed. 



Verma and Denkerl (120111 ) also averaged higher- cadence flow maps over increasing time pe- 
riods — 1, 2, 4, 8, and 16 hr — and found statistics of averaged-fiow speeds (mean, median, 
10th percentile, maximum) all systematically decreased. Such averaging is not the same as 
tracking over longer At, as we have done, but the trends are qualitatively similar: a log-log 
plot of the averaging peri ods and mean speed from tracking granulation given in Table 1 of 
Verma and Denkeii (120111 ) also shows a power-law decrease of mean speed with time between 
1 - 8 hr, though a linear fit to the logarithms returns an exponent near -0.18, about half 
our result. It is plausible that short-lived flows, resolved with shorter At, would primar- 
ily reflect the essentially random effects of granulation, and therefore average incoherently 
over long times. However , comparison of our large-At flow speeds with speeds inferred by 
Verma and Denkerl (120111 ) with a shorter tracking interval At' and then averaged over longer 
time intervals similar to ours implies that shorter-timescale flows are averaging coherently. 
It is possible that the motion tracers — either the intensity pattern, or magnetic flux — used 
in one approach or th e other (or even both !) do not faithfully represent the true plasma 
velocity. In this vein, iRieutord et al.l (120011 ) has argued that the motion of granules does 
not accurately reflect underlying velocities on spatial and temporal scales that are too short. 
This suggests flows inferred from magnetic tracking by flow maps derived with short tracking 
intervals but averaged over long At and flow maps derived with a tracking interval of At 
should be compared in a future study. 

We now return to the question of how flow lifetime r varies with spatial scale. As above, we 
define the spatial scale of an estimated flow as the product of macropixel size and windowing 
parameter a, and therefore have multiple sequences of flow estimates at a few spatial reso- 
lutions for each cadence At. For each sequence of flow estimates, we can estimate the flow 
lifetime r over the whole active region from a one-parameter fit to the autocorrelation as 
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function of lag time AT, assuming exponential decay, exp {—AT/t). In Figure [151 "we plot 
flow lifetimes estimated this way, as a function of spatial scale, for several tracking cadences 
At. In the shorter-scale region of each log-log plot, flow lifetime seems to scale, roughly, in a 
straight line with spatial scale. To characterize this dependence, we fit power-law exponents 
for each cadence over the range of seemingly straight-line log-log scaling (shown with a solid 
line on each plot). For longer cadences, flow lifetime appears to scale nearly linearly with 
scale length, while the dependence of flow lifetime on length scale seems weaker for shorter 
cadences. At larger length scales in each plot, the relationship betwee n lifetime and length 
scale becomes less coherent as binning and apodization are varied. IWelsch et al.l (120091 ) 
found lifetimes at ~ 12 Mm to be two to three steps at their At = 96 min magnetogram 
cadence. For similar spatial scale and At (64 and 128 min) we find shorter lifetimes, by 
a factor of two or so, although this is also the spatial scale where our lifetime estimates 
become incoherent. For all spatial scales we investigate, the flow lifetime is shorter than 
that predicted from a linear scaling of lifetime with length scale, sometimes dramatically so. 
We note that our estimates of lifetimes are rough approximations. Longer-lived flows (e.g., 
supergranular flows) would remain undetected in our lifetime estimates if their asymptotic 
autocorrelation lies below our 1/e cutoff. The longer lifetimes of larger-scale flows is another 
possible explanation of the correlation between flow lifetime and magnetic field strengths: 
stronger magnetic fields might couple nearby regions of plasma more strongly (making the 
plasma more rigid), leading to larger spatial scales for flows in more strongly magnetized 
regions; and such flows should persist for longer than smaller-scale flows. 

In Figure dni we show fitted flow lifetimes for x and y components of the flow (denoted 
+ and X, respectively) as functions of (f), the flow speed averaged over space and time 
for each combination of tracking interval At, spatial binning, and apodization parameter a. 
Plot symbols are color-coded by spatial scale, defined as the product of spatial binning and 
apodization parameter. Generally, higher average speeds correspond to shorter lifetimes. 
While a range of lifetimes was found at each average speed, there appears to be a rough 
upper limit at a given speed, with the peak fitted lifetime scaling roughly as {v)~^. 
Noting that the product of a timescale and the square of speed yields a quantity with the 
dimensions of kinematic viscosity (length^ /time), the mean and standard deviation of r^.iv)'^ 
are (2.7 ± 1.9) x 10^^ cm^ sec~^ and (2.9 ±2.3) x 10^^ cm^ sec~^ for the x and y components 
of the flow, respectively. The limit of flow lifetimes shown by the dashed line in Figure [161 
1.25(f)~^, would correspond to a kinematic viscosity of 7.5 x 10^^ cm^ sec~^. While the 
physical significance of this value is unclear (if it is even physical — it might be an artifact 
of our methods), i t roughly agrees with values of viscosity near the photosphere reported by 



Nesis et al.l fll990l ). 
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5.3. Lifetimes of Curls &: Divergences 

The lifetimes of normal curls and horizontal divergences of the flow flelds, deflned as 

n ■ (V;, X u) (4) 
(V. ■ u) , (5) 

respectivel y, are also of interest: vortica l flows can efficiently inject magnetic energy and he- 



licity (e.g., iLongcope and WelschI 120001). and di verging flows are as sociated with flux e mer- 



gence in both simulations (lAbbett et al.l l2000l ) and observations (IWelsch et al.ll2009l ). In 



Figure [T3| we have overplotted mean absolute values of curls (divergences), averaged over 
tracked pixels in x — y — t datacubes of our flow maps, in red (blue), computed using flnite 
difference approximations with the smallest binning, twice the native 0.16 "pixel size, in the 
denominator. For consistency with our definition of spatial scale used above, the denomi- 
nator should be the product of macropixel size and apodization parameter a. This would, 
however, introduce an additional inverse scaling with the variable plotted on the horizontal 
axis. Average values of curls and divergences at our smallest length sca les for short At are 



on th e order of a few times 10 sec , similar to values reported by IVerma and Denkei 



(I2OIII ). For fixed At, curls and divergences fall by roughly two decades over the approxi- 
mately two decades in length scales we study, implying curl and divergence magnitudes scale 
approximately as 1 /length. Average magnitudes of curls and divergences decreased by about 
a decade over the range of At's in our study, a faster decrease with At than for speeds. 

In Figure [T71 we show autocorrelations of normal curls (red) and horizontal divergences 
(blue) computed in subregions of the full FOV, superimposed on 50 G and 200 G contours 
of I-BlosI) for fiows estimated from 2^-binned data, with cr = 4 macropixels and At = 64 
min. In Figure [TSl we show the fitted lifetimes for curls (left) and divergences (right) ver- 
sus subregion- and boxcar-averaged | -Blob I from NFI, and subregion-averaged \Bz\ and |B| 
from SP, where, as above, lifetimes were determined by one-parameter fits to subregions' 
autocorrelation functions assuming exponential decay. As with the lifetimes of x and y com- 
ponents of the fiow, curl and divergence lifetimes exhibit a dependence on both field strength 
and the median occupancy of above-threshold pixels in each 32^ subregion of 2 x 2-binned 
full-resolution pixels. The rank-order correlation of field strength |B| with curl (divergence) 
lifetime, however, is a statistically significant 0.59 (0.58) for the 59 subregions with occupan- 
cies of 1000 or more (of a possible 1024), implying the lifetime versus field-strength trend is 
independent of any occupancy-lifetime correlation. We repeat that our estimates of |-Blos| 
suffer from non-monotonicity in the Stokes V/I ratio, while our spatial averages of \Bz\ and 
|B| ignore the likelihood of small filling factors for magnetic field structures. Despite these 
limitations, the existence of a trend seems clear, and we find a stronger dependence of the 
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lifetimes of curls on field strength than for lifetimes of divergences. We note that, as with 
the lifetime of flows themselves, the fitted dependence of lifetime of curls and divergences on 
|B| appears weaker than that on either \Bz\ or | -Blob I- 

Finally, in Figure [191 we show the lifetimes of normal curls (red x 's) and horizontal diver- 
gences (blue -|-'s) versus spatial scale, defined as above as the product of macropixel bin 
size and windowing parameter a, for each tracking interval At. Because we estimated flows 
for several macropixel sizes and a's, we have multiple samples at some resolutions, most of 
which roughly agree. At most cadences, the lifetimes increase with spatial scale, and then 
saturate. In Figure [19], we also show linear fits to the logarithms of lifetimes as functions of 
logarithms of spatial scales over the range of spatial scales before saturation occurs. 

6. Summary &: Discussion 

We used FLCT to estimate flows from a ~ 13 hour sequence of high-cadence (~ 2 min), 
high-resolution (0.16 "pixel, 0.3 "diffraction limit) line-of-sight magnetograms of NOAA AR 
10930. This seeing-free, high-resolution, long-duration dataset enabled us to conduct a 
detailed study of the effects of tracking parameter selection on the estimated flows, and to 
investigate the lifetimes of flows and flow structure (normal curls and horizontal divergences) 
as functions of spatial scale and magnetic field strength. For context, we also characterized 
the lifetimes of magnetic structures. 

Several conclusions can be drawn from our results, which may be divided into practical 
conclusions relevant to the process of flow estimation from magnetograms, versus physical 
conclusions about surface magnetic flows and evolution. Our practical conclusions are: 

• Flows estimated by magnetic tracking (and related optical flow) methods with a par- 
ticular set of tracking parameters (i.e., tracking interval At between images, pixel size, 
apodization window size a) are primarily sensitive to evolution on a particular spatial 
scale and time scale. 

• If the tracking interval At between pairs of magnetograms is too short, then changes 
in the field between magnetograms is mostly due to noise instead of actual magnetic 
evolution. This is the noise-dominated regime. Its hallmark is rapid fluctuations in 
estimated velocities, which can be quantified by autocorrelating maps of the velocity 
components. The effect of noise can be reduced by: increasing At, averaging successive 
magnetograms prior to tracking, or both. 

• Because magnetic features persist for much longer than the intensity features to which 
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track ing methods have also been apphed (e.g., iBerger et aLlll998l : IVerma and Denkei 



2OIII ). the upper hmit for At is not strongly constrained in magnetic tracking. If, 
however, the tracking interval At between pairs of magnetograms is longer than the 
flow lifetime at that spatial scale, then the estimated flows reflect an average velocity 
over more than one flow lifetime. This is the displacement-dominated regime. Although 
this time-averaged velocity is physically meaningful, it does not represent the actual 
plasma velocity field at a particular time (e.g., the center of the tracking interval). 

• Hence, choice of At can determine whether flow estimates are noise-dominated (At too 
short), displacement-dominated (At longer than the flow lifetime), or representative 
of actual flows. To estimate actual plasma velocities, appropriate choice of At can be 
established through post-facto analysis of flow maps, by autocorrelating flows to find 
a At that both maximizes frame-to-frame autocorrelations (which thus minimizes the 
effects of noise) and is less than the flow lifetime r, which we take to be the time at 
which the autocorrelation drops to 1/e. 

• Average estimated speeds decrease with increasing tracking interval At between mag- 
netograms (see Figure [li" 

Our physical conclusions are: 



As noted by previous researchers (e.g.. iBerger et al.lll998[ ). photospheric flows operate 



over a range of spatial scales, and exist for a range of lifetimes. Hence, referring to 
"the flow field," by itself, is imprecise: explicit reference to spatial and temporal scales 
must also be made. 

Magnetic structure persists substantially longer in stronger-field regions generally, and 
longer in regions of strong vertical field \Bz\ than in regions of strong magnetic field 
strength |B|. 

Magnetic structure persists substantially longer when averaged over larger spatial 
scales. Flow components and their derivatives estimated on larger spatial scales also 
persist for significantly longer, though lifetimes scale less than linearly with length 
scale. 

Flow components, horizontal divergences, and normal curls persist somewhat longer in 
regions with stronger fields: our estimates of their lifetimes are significantly statistically 
correlated with spatially averaged values of each of |5los|, \Bz\, and |B|. We find the 
lifetimes of curls exhibit a slightly stronger dependence on field strength than either 
flows components or horizontal divergences. 
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• Flows with faster average speeds typically exhibit shorter lifetimes. Also, there appears 
to be a rough upper limit on flow lifetime at a given average speed, which scales as the 
inverse of the square of the average speed. 

In this study, we found sets of differing velocity fields for a single magnetogram dataset 
by applying a single correlation tracking method over various spatial and temporal scales. 
This is not surprising: coherent flow structures on varying scales are known to be present 
in the photospheric flow field, e.g., global-scale Sun's differential rotation (in our synodic 
frame) and meridional flow, high-speed granular flows and lower-speed supergranular flows 
driven by convection, and possibly flows arising from Lorentz forces in magnetized regions. 
In addition, turbulent flows are likely present. This, however, raises the question, "Given 
that a family of velocity fields can be derived on several scales, how can one determine the 
best estimate of the velocity field?" But what determines "best?" It is tempting to try 
to estimate the "true" velocity field in each resolution element (pixel), but what does this 
mean physically? Assuming the photospheric plasma behaves exactly as a fluid down to 
some kinetic scale, then within our smallest resolution element there is a distribution of 
subresolution velocities, and a meaningful definition of the "true" velocity must be defined 
in terms of this distribution. If the distribution of velocities is peaked in speed and direction, 
then some statistical property of each (e.g., mean or median) might meaningfully quantify 
the flow. It is, however, probable that sub-resolution flows are consistent with a power-law 
distribution of flow energy with spatial wavenumber, similar to that proposed by Kolmogorov 
(though perhaps with a different wavenumber scaling), through some inertial range down 
to a dissipation s cale. A reasonable expectation based upon o ur results, as well as those of 



previous studies (iBerger et al.l Il998l : IVerma and Denkerl l201ll ) , is that flows on successively 



smaller scales will have higher speeds. In this case, what we infer from the motion of 
magnetic flux (or of brightness variations, if tracking intensity) is an "effective flow" — the 
combined effect of flows acting on all scales. But it is probable that the effective flow does 
not correspond to "the" physical velocity. 

From the practical standpoint of understanding how flows are related to other phenomena 
of interest (e.g., flares, CMEs, coronal heating, the dispersal of active region magnetic flux), 
we must investigate relationships between these phenomena and flows on varying spatial and 
temporal scales to understand which (if any) scales are most relevant. A key challenge is that 
the largest velocities are present at the smallest scales: while in a sense these strongest flows 
are the easiest to detect, they also have the shortest lifetimes, and are therefore probably the 
least relevant for some processes with longer characteristic time scales, e.g., the emergence 
of large-length-scale magnetic flux systems in active regions. Computing the Poynting flux, 
for instance, based upon an instantaneous estimate of flows on granular scales will likely 
give an incomplete (if not misleading) picture of how much and where magnetic energy is 
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being transported across the photospheric layer. Hence, much work remains to understand 
photospheric flows and their relationships to processes like the transport of magnetic energy 
across the photosphere. 

We remark that the techniques we have applied here can also be applied to both real and 
simulated vector magnetogram sequences, to better understand evolution of both magnetic 
fields and flows. Long sequences of vector magnetograms (from, e.g., HMI) could be tracked 
to analyze flow properties over a range of spatial and temporal scales as functions of the 
vector magnetic field's properti es. The lifetinaes of magnetic and flow fields from simulations 
of photospheric evolution (e.g., ICheung et al.ll2010l ) can also by analyzed via autocorrelation 
for comparison with lifetimes inferred from observations. 

We e xpect observ a tions using the ne west generation of gr ound-bas ed solar telescop es 



NST Goode et al. 



J2010h. GREGOR dVolkmer et al.lboioh . ATST faimmele et"aDl201oh . 
and EST (IZuccarello and Zuccarelldl201ll ) — to play a key role in future studies of flows at 



high resolution. 
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A. Calibration by Co-alignment with SP Vector Magnetograms 



Schrijver et al.l (120081 ) used SP data to produce vector magnetograms of AR 10930 over 12 - 
13 Dec. 2006, which are available online. One of the magnetograms they produced overlaps 
with our tracking run. Here, we describe the procedure we used to interpolate magnetic field 
data from that magnetogram onto a downsampled, cropped grid corresponding to our NFI 
observations. Note that our ultimate goal is to be able to relate the lifetimes of magnetic field 
structures and flows over regions of the NFI FOV to average magnetic properties in those 
regions. Hence, we do not seek pixel-scale agreement between the SP and NFI magnetograms. 
Rather, we seek approximate agreement on the scales of a few arcseconds. 
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Schrijver et al.l (120081 ) report sampling the SP data into 0.63 "pixels. We resampled NFI 
observations of -Blos at the start, middle, and end of the 45-minute SP scanning interval from 
the original 0.16 "pixels into (4 x 4) bins to achieve commensurate resolution. Comparison of 
features in the Bz map from SP (Figure [201 upper left) with binned -Blos maps from NFI at 
the start and end of the scan (Figure |20l middle and bottom of left column) shows that most 
of the discrepancies between the SP and NFI data are in the horizontal direction: features 
in the SP image extend over N^^sp = 295 pixels in the horizontal direction, while the NFI 
images of the corresponding features are only Nx^nfi = 240 pixels in horizontal extent. 

Given the position of the AR near 20° W from disk center at the time of the scan, and the 
~ 10° extent of the AR, the cumulative effect of horizontal distortion in pixel lengths from 
foreshortening in NFI pixels compared to reprojected SP pixels over the FOV should account 
for a ~ 9% difference, or about 22 pixels. Hence, additional distortion is present. We also 
note that, even if reprojection fully accounted for the SP/NFI discrepancies, un-reprojecting 
the SP data would require a difficult inversion of the convolved effects of solar rotation, 
spacecraft pointing, and spectrograph rastering over the 45 minute scanning period. 

We therefore adopt an ad hoc approach, using repeated Fourier interpolations to spatially 
interpolate the SP data onto each column of the NFI grid. We first interpolated the SP 
data onto N^^nfi evenly-spaced shifts Nx,sp/Nx,nfi apart (one Fourier interpolation per 
shift), but found systematic residual horizontal distortions were present, consistent with a 
higher-order variation in distortion across the image. We therefore included a quadratic 
term in the spacing between interpolations: where is the location of the i-th column of 
NFI data, we interpolated the SP data at x[ = Nx,sp/Nx,NFi{,o,Xi + 6xf ). By trial and error, 
we found a = 0.925 and b = 0.06 gave good qualitative agreement between ±100 G and 
±400 G contours of -Blos from NFI at the middle of the scan interval with from the SP 
magnetogram, as shown in the upper right panel of Figure I2U1 It may be seen that the NFI's 
contours of -Blos do not match SP's B^ everywhere, but the remaining offsets do not appear 
to be systematic across the image, or along columns. (We note that the non-monotonicity in 
the NFI signal precludes use of some quantitative measures of NFI-SP agreement [e.g., cross- 
correlation, or minimized square difference] since optimal alignment probably corresponds 
to the weakest umbral field in NFI being co-aligned with the strongest umbral field in SP, 
while features outside of the umbra should generally correspond one-to-one.) 

As discussed above, we are primarily interested in average magnetic properties (vertical 
field strength, total field strength) in larger-scale regions of the NFI FOV, so we only use the 
interpolated SP data binned into (8 x 8) and (16 x 16) macropixels. Comparisons of (32 x 32)- 
binned NFI |-Blos| data with (8 x 8)-binned |-B^| (and |B|) from SP in the middle-right (and 
bottom-right) panel of Figure [20] shows good spatial agreement at that resolution. 
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Note that the SP data do not fully cover the NFI FOV, which is 256^ pixels when binned 
(4x4). Consequently, we can only use the SP data on a subset of the full NFI FOV. 
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Fig. 4. — The grayscale image is |i?LOs| in a full-resolution magnetogram, with saturation at ±4 kG. Red 
curves show rank-order autocorrelations of i?LOS, in (32 x 32)-binned subregions of (2x2)- binned data, 
vs. lag time, up to the duration of the dataset. The vertical range in each cell is [-0.5, 1.0], with a dashed 
black line at zero correlation. Blue curves show one-parameter fits to the decorrelation in each subregion, 
assuming exponential decay, with the decay constant as the only free parameter. Only subregions with 
median occupancy of at least 20% of pixels above our 15 G threshold were Gt. As expected, held structures 
persist longer in stronger-held regions. 
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Fig. 5. — Fitted lifetimes of boxcar-averaged |i?LOs| s-nd |B|, averaged over each subregion, for the non- 
bubble region of the FOV. Blue triangles correspond to (32 x 32)-binned subregions (from fits in Figure\^ of 
I^LOS I • Rod + 's correspond to (16 x 16)-binned subregions |i?LOS I (autocorrelations and fits at this resolution 
are not plotted in this paper), and green x 's correspond to |B| averaged over (16 x 16)2-binned subregions. 
Fitted exponents to log-log values of ji^LOsI and |B| for (16 x 16) bins are shown with corresponding colors. 
Dashed lines correspond to uncertainties in each fitted slope, based upon x^; the slope of the blue fit (0.68) 
is similar to the red fit (0.69). 
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Fig. 6. — Plotted lines show frame-to-frame correlation coefficients for the x- and y- components of estimated 
Bows, for varying time intervals At between tracked magnetograms — 2 min in black, 4 min in aqua, 8 min 
in blue, 16 min in red. Dashed curves are from Rows derived from un-averaged magnetograms, while solid 
curves are from Rows derived from boxcar-averaged magnetograms. Flow estimates are more stable with 
averaged magnetograms and with longer At, until saturation occurs around At ^ 8 min. 
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Fig. 7. — Successive Row maps in a subregion of the full field of view, derived from unaveraged magnetograms 
(binned 2x2) separated in time by At — 2 min with a = 4 pixels. Initial flows are shown by red vectors, 
subsequent flows are shown by blue vectors. Note significant variations between the flow fields in this 
noise-dominated regime. 
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Fig. 8. — Successive Row maps from the same subregion, derived from averaged magnetograms (binned 2x2) 
separated in time by At — 8 min with (7 = 4 pixels. Initial Hows are shown by red vectors, subsequent Bows 
are shown by blue vectors. Note signihcant agreement between Bow Gelds. 
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Fig. 9. — Plotted lines show correlation coefRcients for the x- and y- components of flows estimated at ti 
and ti + AT, for several time intervals At between tracked magnetograms — 4 min in black, 8 min in aqua, 
16 min in blue, 64 min in red. Dashed curves are from flows derived with a — 4 pixels, while solid curves 
are from Rows with (7 = 8 pixels. All flows were derived from boxcar-averaged magnetograms. Horizontal 
dotted line is 1/e. 
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Fig. 10. — Top: Autocorrelations of x- and y-components of flows derived using two values of a (dashed = 
4 pix., solid = 8 pix.) from magnetograms binned into different-sized macropixels (black = 2^, Wue = 4^, 
red = 8^). Evidently, keeping a fixed but doubling bin size approximates keeping bin size fixed but doubling 
(7. Additional binning did not improve correlations. Bottom: Autocorrelations of x- and y-components 
of flows estimated with large values of a (16 pix. for all flows here) and from data rebinned into larger 
macropixel sizes (colors as in top panel) for At = 128 min (solid) and At = 256 min (dashed) are long-lived. 
As in top panel, additional binning did not improve correlations. All flows in both panels were derived from 
boxcar-averaged magnetograms. 
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Fig. 11. — The grayscale image is IBlosI in a full-resolution magnetogram, with saturation at ±4 kG. 
Autocorrelations of Ux (red) and Uy (blue) in subregions of the full field of view. The flows were derived 
with At = 64 min and a = A pix with (2 x 2) binned data. The vertical range in each cell is [-0.5, 1.0], 
with a dashed black line at zero correlation. Autocorrelations are essentially random in weak- field regions, 
in which few pixels were above the threshold for tracking. 
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Fig. 12. — Fitted How lifetimes for non-bubble subregions of Figure [711 in which at least 20% of the 32^ 
(2 X 2) pixels exceeded the tracking threshold, as a function of the spatial average of the unsigned, boxcar- 
averaged LOS magnetic Geld in each subregion. Lifetimes for Ux in subregions are shown with red -h's, Uy 
are blue x 's, and Gts to these are red and blue solid lines. The red A's (green Ws) are lifetimes of 
versus subregion-averaged \Bz\ (\'B\) from the SP data, and the black (green) line is a fit, with slope 0.17 
(0.13). The correlation between flow lifetime and field strength is statistically significant, though the fitted 
power-law dependence is weak. 
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Fig. 13. — Black numerals show average How speed (left axis) vs. spatial scale, parametrized by the product of 
the macropixel bin size and the windowing parameter a, for several flow cadences; the numbers correspond 
to log2 (bin size), in units of 0.16 "pixels. Red (blue) numerals show the log of the average normal curl 
(horizontal divergence) of the flow, on the right scale, and are slightly offset right (left) for clarity. Because 
we estimated flows for several macropixel sizes and a's, we have multiple samples at most resolutions. For 
shorter At, bin size affects the average How speed at a Hxed spatial scale; for longer At, bin size does not 
strongly affect average How speed at a Hxed spatial scale. In contrast, binning seems to affect average How 
curls and divergences more at larger At. 
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Fig. 14. — Flow speed as a function of cadence, averaged over both space and time for each spatial resolution 
(binning xa) and over spatial resolutions at each At. Error bars show standard deviation over resolutions. 
Least squares linear fit to the logarithms of cadence and speeds returns a power-law exponent of -0.34. 
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Fig. 15. — Flow lifetimes (+'s for Ux, x's for Uy) as a function of flow spatial scale, as parametrized by 
the product of the macropixcl bin size and the windowing parameter a , for several Row cadences. Because 
we estimated Rows for several macropixel sizes and a 's, we have multiple samples at some resolutions. We 
computed linear fits (solid lines) to the logarithms of spatial scale and lifetime over a limited range of scales. 
The fitted power-law exponent is printed in the title of each cadence's plot. 
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Fig. 16. — Fitted lifetime vs. How speed averaged over space and time for each choice of At, spatial binning, 
and a. Lifetimes for and Uy are plotted with + 's and x 's, respectively. Generally, higher average speeds 
correspond to shorter lifetimes. While a range of lifetimes exists at each average speed, there appears to be 
an upper limit at a given average speed, with the peak fitted lifetime scaling roughly as (average speed)^^ 
(dashed line; note that this line is not a fit). Points are color-coded by spatial scale of the flow (binning xa). 
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Fig. 17. — The grayscale image is IBlosI in a full-resolution magnetogram, with saturation at ±4 kG. 
Autocorrelations of normal curls (red) and horizontal divergences (blue) in each subregion of the full field 
of view are overplotted. The Rows were derived with At = 64 min and cr = 4 pix. The vertical range in 
each cell is [-0.5, 1.0], with a dashed black line at zero correlation. Autocorrelations are essentially random 
in weak-field regions, in which few pixels were above the threshold for tracking. 
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Fig. 18. — Fitted lifetimes of normal curls (left) and horizontal divergences (right) for subregions of Figure 
[T7\ in which at least 20% of the 32^ (2 x 2) pixels exceeded the tracking threshold, as functions of the 
spatial average |i?LOs| (from NFI; blue +'s), \Bz\ (from SP; red A'sJ, or |B| (from SP; green Vs) in each 
subregion. The correlation between curl/ divergence lifetime and field strength is statistically significant, 
though the fitted power-law dependence is weak. The magnetic dependence of the lifetimes of curls is 
somewhat stronger than that of lifetimes of divergences. 



-47- 



At = 2 min. 



At = 4 min. 



1 10 



> 

b 1 



^iv. fit = 0.13 




: >^ 










+ 




X 




X 




djrl fit = 0.048 



1 10 100 

Spatial scale, binsize * a [Mm] 

At = 8 min. 



1 



Div. fit = 0.40 




X 



Curl fit ^ 0.56 

1 10 100 

Spatial scale, binsize * a [Mm] 

At = 32 min. 



- 1 



Div. fit = 0.55 




X 



Curl fit = 0.66 

1 10 100 

Spatial scale, binsize * a [Mm] 

At = 128 min. 



100 



o 1 



Div. fit = 0.27 



1 



X )t( 



X X 



> 

Q 1 



Curl fit ^ 0.21 

1 10 100 

Spatial scale, binsize * a [Mm] 

At = 1 6 min. 



+ "Div. fit 



10 - 



0.49 




X 
X 
X 



X 



Curl fit = 0.50 

1 10 100 

Spatial scale, binsize * a [Mm] 

At = 64 min. 



^- lOOrOiv. fit = 0.54 



D 10- 




Curl fit ^ 0.73 

1 10 100 

Spatial scale, binsize * a [Mm] 

At = 256 min. 



Div. fit = 


0.57 


X ¥ 


X 

4^ 








e 










: + 




+ + 






-""^ X 
























Curl fit 


= 0.76 : 


> 


1 


10 


100 




Spatia 


scale, binsize * 


a [Mm] 





X 1000 r 



:Div. fit = 0.66 



1 r 




Spatial scale, binsize * a [Mm] 



Fig. 19. — Lifetimes of normal curls (red x 's) and horizontal divergences (blue +'s) as a function of flow 
spatial scale, as parametrized by the product of the macropixel bin size and the windowing parameter a, 
for several flow cadences. Because we estimated Rows for several macropixel sizes and a 's, we have multiple 
samples at some resolutions. We computed linear fits (solid lines) to the logarithms of spatial scale and 
lifetime over a limited range of scales, and list the fitted power-law exponents on each cadence's plot. 
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Cropped SP, 2006 Dec. 12, 20:30 - 21:15 Grayscale: Interpolated SP; Countours: NF 




Fig. 20. — Left column: a) SP map of Bz, scanned over 45 minutes, b) Sz c) NFI images binned to (4 x 4) 
from 0.16 "pixels, from the beginning and end of the SP scan. Note lateral extent of SP data is 295 pixels, 
while NFI data are 240 pixels. Hence, the primary difference between the datasets is lateral stretching. Right 
column: d) Contours of Blos from scan-centered NFI magnetogram (—/+ lOOG and 400 G in red/blue, 
rcsp.) overlain on the interpolated SP data. Slight inaccuracies in co-alignment vary across the FOV. e) &: 
f) |i?LOsl from NFI contoured over grayscale of\Bz\ and |B| from SP, resp., both binned (8 x 8). Contours 
of |Blos| from NFI correspond to 100, 200, 400, 800, and 1600 G. 



